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Motivated by the absence of dark matter signals in direct detection experiments (such as the 
recent LUX and PandaX-II experiments) and the discovery of gravitational waves (GWs) at aLIGO, 
we discuss the possibility to explore a generic classes of scalar dark matter models using the 
complementary searches via phase transition GWs and the future lepton collider signatures. We 
focus on the inert scalar multiplet dark matter models and the mixed inert scalar dark matter 
models, which could undergo a strong first-order phase transitions during the evolution of the early 
universe, and might produce detectable phase transition GW signals at future GW experiments, 
such as eLISA, DECIGO and BBO. We find that the future GW signature, together with circular 
electron-positron collider, could further explore the model’s blind spot parameter region, at which 
the dark matter-Higgs coupling is identically zero, thus avoiding the dark matter spin-independent 


direct detection constraints. 


PACS numbers: 


I. INTRODUCTION 


A brand new door has been opened to study the fun- 
damental physics and the particle cosmology by the grav- 
itational waves (GWs) after the discovery of the GWs 
by Advanced Laser Interferometer Gravitational Wave 
Observatory (aLIGO) [I]. One active research field is 
the idea of probing the new physics through phase tran- 
sition GWs, where a strong first-order phase transition 
(FOPT) is induced by the new physics models and can 
produce detectable GW signals from three mechanisms: 
collisions of expanding bubbles walls, magnetohydrody- 
namic turbulence of bubbles and sound waves in the hot 
plasma of the early universe [2H9]. Various new physics 
models have been studied by the new approach of GW 
signals in the post aLIGO era [I0}27], such as the GW 
detection of the dark matter (DM) [24], the hidden 
sector [11] [12] [9] [22] and the electroweak baryogene- 
sis [O07 S 0 27. 


Another important puzzle in particle cosmology is the 
particle nature of the DM. Usually there are three ways 
to search for particle DM: direct detection, indirect de- 
tection, and collider searches. With the experimental 
precision of the DM direct detection are gradually ap- 
proaching the neutrino backgrounds, such as the recent 
LUX and PandaX-II experiments [28] [29], GWs may be- 
come a new approach to explore the existence of DM 
since in a large classes of DM models where a FOPT 
can be trigged by the DM particles and other associ- 
ated particles [30}32]. In general, a strong FOPT can be 
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triggered by the DM candidates in generic new physics 
models, which can produce detectable GW spectrum by 
Evolved Laser Interferometer Space Antenna (eLISA) [33], 
Big Bang Observer (BBO) [34], Deci-hertz Interferome- 
ter Gravitational wave Observatory (DECIGO) and 
Ultimate-DECIGO (U-DECIGO) [36]. 


In this paper, we study the chance to explore a general 
classes of DM models through phase transition GW signals 
since the GW signal becomes a novel approach to study 
the property of the various new physics models with 
extended Higgs sector after the discovery of Higgs boson at 
LHC and GWs at aLIGO. And the extended Higgs sector 
could incorporate the scalar DM candidate. One generic 
classes of new physics models are the inert multiplet scalar 
DM models [87] [38] and mixed scalar DM models [39] [40]. 
In these models, the scalar multiplets usually belong to a 
hidden sector, and thus contain the DM candidate. At 
the same time, we know that usually extra scalars could 
change the phase transition structure for the standard 
model (SM) Higgs boson around the temperature of the 
electroweak scale. Due to extra scalar degree of freedom in 
the thermal plasma, the extra inert scalars could enhance 
the strength of the electroweak phase transition. On 
the other hand, these models typically encounter the 
tension between tight direct detection constraints and 
the strong FOPT. To avoid the tighter and tighter direct 
detection constraints, we could focus on the blind spot 
region, at which the Higgs-DM coupling is identically zero, 
thus avoiding the DM spin-independent direct detection 
constraints. In the scalar multiplet DM models, usually 
there is strong correlation between the Higgs-DM coupling 
and the strong FOPT: zero Higgs-DM coupling indicates 
there is no effect on the electroweak phase transition from 
the inert DM sector. The only exception in the single 
scalar multiplet DM models is the inert doublet model 


(IDM) [41] [42], in which although the direct detection rate 
could be tiny the strong FOPT could be realized. In the 
blind spot region, requiring the observed relic abundance, 
the allowed parameter space in IDM is very limited [43]. 
To enlarge the allowed parameter region, we consider 
the mixed inert DM models, in which two inert scalars 
are mixed [39] [40], and thus there is large blind spot 
parameter region. We will study the possibility to obtain 
strong FOPT in the blind spot parameter region in the 
mixed DM models. We expect that the GW signature 
could further explore the parameter region which has not 
yet been explored by the direct DM detection experiments. 


This paper is organized as follows: In Section I, we 
discuss the general inert scalar DM model, and point out 
the IDM in the blind spot region could lead to the GW 
signals while avoiding the DM constraints. In Section 
III, the mixed inert scalar models are investigated in the 
blind spot region with the DM constraints and the GW 
signals. In Section IV, we show our final discussions and 
conclusions. 


II. MULTIPLET DARK MATTERS AND 
GRAVITATIONAL WAVES 


Let us start from a class of DM models in which the 
SM is extended by adding an electroweak scalar multiplet. 
The scalar multiplet is Z2-odd under an imposed global Zə 
symmetry. In certain hypercharge assignment, the neutral 
component in the scalar multiplet could be the lightest 
particle in the multiplet, and thus it could be a possible 
DM candidate. This class of models have been well- 
studied in the inert isospin-singlet [44], doublet [42], 
and triplet cases. And a general DM multiplets 
in terms of the SU(2); representations have also been 
investigated systematically [87] [38]. In general, the scalar 
multiplet Hn belongs to the representation n of the SU (2) 
group, with hypercharge Y under the transformation of 
the SU(2), x U(1)y gauge group 


hti 
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where the T? charge j runs from —25+ to 454+. The 
relevant Lagrangian is written as 


L= D, H} D" H, — V (©, Hn), (2) 


where ® is the SM Higgs doublet. The scalar potential is 
usually taken to be 


V(®, Hn) = V(®) + M2H? + doHA +A (Hir? Hn)? 
+ 30H? + ALDIT O) (Air? Hn), (3) 


where V(®) = —p?®? + At. For a real multiplet t, 
the term HÌ r” Hn is identical to zero, and thus the A5 
and à}, terms disappear. As a consequence, only the 
term A3 connects H, to ®, and the tree-level masses 
of the multiplet components are degenerate with m2 = 
M2+ Agu? 2. Similarly, in the complex singlet model, the 
à% and à% terms do not exist. For complex doublet model, 
there will be more terms than Eq. in the potential 


V(%, H2) = V(®) + MZH? + `H 
+ à? H? + Að H|? 
+ 5/2 [(®'H2)? + h.c.] . (4) 


To make the real component H as the DM candidate, 
A345 = Ag + A4 + à5 < Az and às < 0 is needed. For a 
complex multiplet with higher isospin, for simplicity we 
will only focus on the case with hypercharge Y = 0 ? in 
this work. As shown in Ref. [38], typically the complex 
triplet with Y = 0 can be viewed as two real triplets. 

Firstly, we discuss the FOPT induced by the inert scalar 
models. Up to one-loop level, the effective potential at 
the finite temperature can be written as 


Vow (®, T) = Vo(®) a Vow (®) oF Viner(®, T) + Vaaisy (®, T), 


where Viher(®, T) + Vaaisy(®, T) is the thermal effects 
with the daisy resummation and Vow(®) is the Coleman- 
Weinberg potential at zero temperature. In the inert 
scalar models, the necessary potential barrier for the 
strong FOPT in the effective potential at the finite tem- 
perature origins from thermal loop effects, where the 
bosons contribute to the effective potential with the form 
Vp > (—T/12m)(m2oco(®,T))?” in the limit of high- 
temperature expansion. For the heavy fields whose masses 
are much larger than the critical temperature, the con- 
tributions from heavy particles can be omitted due to 
Boltzmann suppression. This can help to simplify our 
discussions when the models have many new fields at 
different energy scales. Thus, due to the above consider- 
ations and the fact that we only study the inert scalar 
models and the thermal barrier induced strong FOPT, 
the parameter space of very large mass scalar new bosons 
are not favored since the field independent term in the 
thermal masses should be small enough in order not to 
dilute the cubic terms. We focus our study in the case 
without very heavy scalars. Therefore, to begin with the 
concrete prediction of the GW signals in the following 


1 A real multiplet only has integer isospin. Thus there is no real 
doublet, quadruplet, etc. 

2 The loop corrections to the tree-level masses cause small mass 
splitting between charged components and neutral DM candidate. 
If the hypercharge Y is not zero, the neutral components of 
the complex multiplet have interaction with the Z boson, which 
usually leads to a large spin-independent cross-section between 
DM and the nucleon {38}. But if for some reason there is a mass 
splitting between real and imaginary components of the neutral 
scalar, the direct detection constraints might be avoided. 
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examples, the general effective potential near the phase 
transition temperature can be further approximated by 


T A 
Peon, (5) 
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Here, h represents the Higgs boson field in the unitary 
gauge as (®) = h/v2, where the angle bracket means 
the vacuum expectation value (VEV) of the field. The 
coefficient € quantifies the total interactions between the 
light bosons and the Higgs boson. In the SM, we have 
ESM = Omie a and csM = Doe 

Thus, in this case of qualitative analysis, the washout 
parameter can be obtained as §v) (Fe) ~ gy, and it should 
be larger than one for a strong FOPT. In inert scalar 
models, the washout parameter usually is approximately 
proportional to the Higgs-inert scalar coupling 3 for the 
potential in Eq. (3). Roughly, a larger Higgs-inert scalar 
coupling produces a stronger FOPT. For high multiplets, 
the situation becomes more complicated due to the large 
thermal mass coming from gauge interactions. The large 
thermal masses on the new scalars cause plasma screening, 
and thus decrease the strength of the phase transition. 
Therefore, there are two parameters which control the 
strong FOPT: the Higgs-inert scalar coupling, and the 
thermal mass coming from gauge interactions. 

We know that the Higgs-inert scalar coupling is usually 
constrained by the DM direct detection experiments. The 
tree-level interaction between the DM and the nucleon 
is through the Higgs boson exchange, which induces the 
elastic spin-independent cross section 


2 m2 2 
ost = fx, 2 (2) ; (6) 
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where the Higgs-DM coupling Any, = A345/2 for inert 
doublet, and Any. = A3/2 for other intert multiplet. Cur- 
rent DM direct detection experiments constrain 
that the Any, is around 0.012 for about 100 GeV DM 
mass. In general, the coupling Any, in direct detection 
is the same as the one which mainly controls the strong 
FOPT. Thus the tighter direct detection constraint, the 
less significance of the FOPT. The only exception is the 
IDM. In the IDM, due to the extra interaction terms in 
the scalar potential, the strong correlation between direct 
detection constraints and FOPT can be avoided. We will 
focus on the blind spot region in IDM, in which although 
the direct detection constraint can be evaded, the strong 
FOPT can be induced. Thus, the corresponding phase 
transition GWs can be produced and the DM abundance 
can be explained. 


A. Inert Singlet, Triplet and Multiplet Scalar 
Models 


In the inert singlet model, the Higgs portal term 
A3|®|7,S2/2 is crucial, where we use S represent the inert 


singlet scalar instead of Hı. The corresponding washout 
parameter in the inert singlet model is approximately 


3 3 
(v)(Pe) SPE + (5/2)°? 
T, 6TÀ Í 


(7) 


In the inert triplet model, we only consider a simple model 
with an SU(2), triplet scalar H3(1,3,0) with a zero hy- 
percharge. The relevant term involving the triplet scalar 
Hz should be A3®'®Tr(H?). The corresponding washout 
parameter in the inert triplet model is approximately 


(oTa) _ ERTE 4 3 (03/2)? 
T, a 6TÀ : (8) 


However, in both models the direct DM search puts strong 
constraint A3 < 0.01, which means the FOPT is forbidden 
by the DM direct experiments. 

For a general n-multiplet inert scalar Hp models, 


there are two competing sources which affect the strong 
FOPT [46]: 


e Higgs-inert multiplet coupling with enw ~ 
n(As/2)"”?; 


e plasma screening due to large thermal mass coming 
from gauge interactions. 


Typically the higher multiplet, the larger screening and 
decoupling effects, which weaken the FOPT. According 
to Ref. [46], for a multiplet with n > 3, the screening 
effects significantly decrease the strength of the FOPT. 
Furthermore, another severe constraint for high multiplet 
model is from Higgs diphoton rate with all the charged 
scalars running in the loop. For the real scalar multiplet 
with n > 3, the Higgs coupling measurement data put very 
strong constraints on the masses of the scalar multiplet to 
be greater than 300 GeV, which makes the scalar degree 
of freedom decoupled from the plasma. 


B. Inert Doublet Model 


As mentioned above, usually the blind spot region with 
zero Higgs-DM coupling indicates the DM sector does 
not affect the electroweak phase transition. In IDM, zero 
Higgs-DM coupling does not indicate zero Higgs-inert 
scalar couplings are zero. There is no correlation between 
direct detection and strong FOPT in IDM. Therefore, we 
expect to obtain the strong FOPT and detectable GW 
signals in the blind spot parameter region. 

We investigate the finite temperature effective poten- 
tial and discuss conditions of strong FOPT in detail 
[49]. The relevant scalar potential in the IDM is given 
in Eq. (4) where Hə stands for the inert doublet scalar 
without VEV. In the IDM, we assume that only ® can ac- 
quire VEV, namely T = (0,v + h)/V2 and H is the 
lightest component of the inert doublet Hə with the 
mass mz, = Mj, + $A345v?. Thus, the particle H is 


the DM candidate here. The other neutral scalar mass is 
m? = MÈ + 3(A3 + A4 — A5)v?, and the charged scalar 
mass is m? = Mĝ + $A3v?. The thermal phase transi- 
tion with full 2-loop effective potential has been studied 
recently and it shows the one-loop effective potential 
in the high temperature expansion is rather reliable in 
the IDM. To clearly see the phase transition physics and 
simplify the following discussions on the phase transition 
GW signals, we take the following approximation of the 
one-loop effective potential including the daisy resumma- 
tion: 


1 AÀ 
Velh, T) ~ 3 (—W? + crpm T”) h? + ra 


T 
— -Ene(mp(h,T))?? 
127 
mp (h T) 5 m?(h,T) 
6472 Ca T? 
mf(h) meh) 


np 


where log cg = 5.408 and log cp = 2.635. In the effective 
potential, the particles running in the loop are the parti- 
cles in the model with the following degrees of freedom: 


nwt 4, nz 2, Nr 3, 
—12. 


Nh = NH = Ny Ny 1, n 


The field-dependent masses of the gauge bosons and the 
top quark at zero temperature are given by 


24 7/2 2 
g +g y 
= 3 h?, m? (h) = 3, 
where y; is the top Yukawa coupling. The field-dependent 
thermal masses at the temperature T are 


2 g9? 2 2 
mẹ(h) = T” , mz(h) 


mi(h,T) = m2 x 3A =p? + aT’, 


m2,(h,T) 5s + dg + Ag)h? + M3, + cT °, 
ETIR 50s + dg — Ag)h? + M3, + cT °, 
m? (h, T) x Saah? +M? toT’, 
where cy = A 2As tA | Sota" | ue and c2 = a2 + 


2-32 
2AstAa 4 3g a. In the above formulae, we have con- 


12 
sidered the contribution from daisy resummation, which 
reads as 


V> oe ([m2(n, T]? — [m2(ny]®”) 


Here, the thermal field-dependent masses m?(h,T) = 
m?(h) + Ui (h,T), where II;(h, T) is the bosonic field 7’s 
self-energy in the IR limit. As in the SM, this cubic term 
is the unique source to produce a thermal barrier in the 
effective potential, and in the Higgs sector extended mod- 
els, the new degree of freedoms in the inert scalar models 


increase the barrier and hence produce strong FOPT. How- 
ever, the cubic terms should be large enough to produce 
a strong FOPT. To avoid diluting the cubic contribution 
to the thermal barrier, the Higgs boson field independent 
term needs to be very small [31)[51]. Thus, in this limit, we 
A3) 3/2 Aa+à4— às ) 3/2 Asas \3/2 
Gey e Oe) 


have EIDM © ESM + 2 
and cDM œ% csm + Then, the corresponding 
washout parameter in the IDM is roughly 


(v)(Te) em tA + aa 4 (Agua) 9/? 


Te 6TA 


From Eq. (6) and the washout parameter, we find that 
small DM direct detection rate and strong FOPT can 
be realized when we take the blind spot region, in which 
Anxx = 345/2 approaches to zero 4. The DM model of 
IDM needs to satisfy the required DM relic density ob- 
served from Planck: Qh? = 0.1184 + 0.0012 [52]. This put 
very strict constraints on the IDM: the DM mass is deter- 
mined to be my > 540 GeV [38] except for the parameter 
region with large mass splitting between charged and 
neutral components. For the DM mass lower than m,,/2, 
the Higgs invisible decay puts very tight constraint on 
the parameter space. According to the latest study [43], 
there are two viable mass regions: 


e near Higgs funnel region with large mass splitting be- 
tween charged and neutral components: m, around 
55 ~ 75 GeV with A345 < 0.04; 


e heavy DM region: my > 540 GeV with A345 in a 
broader range as m, gets heavier. 


To keep the scalar non-decoupled from the thermal plasma, 
it is necessary to have light DM. The dominant DM anni- 
hilation channel will be yx > WW*, ZZ* with contact, 
t- and s-channels. We will focus on the DM mass around 
55 ~ 75 GeV, and the blind spot region with A345 ~ 0. 
Combined the direct DM constraints, the DM relic den- 
sity, collider constraints [43] and the conditions for strong 
FOPT, this light mass region 55 ~ 75 GeV is favored. The 
strong FOPT can be produced if A3/2 and (A3+A4—A5)/2 
are order 1, then detectable GW signals can be produced, 
while keep the coupling between Higgs boson and DM 
pair small enough to satisfy DM direct experiments and 
relic density. 

Considering the above discussion, we take one set of 
benchmark points A3 = 2.84726, A4 = A5 = —1.41293 
and Mp = 60.89 GeV. Then, the corresponding DM 
mass is 66 GeV, the pseudo scalar mass and the charged 
scalar mass are both 300 GeV, Anyy = A345/2 = 0.0107. 
Taking this set of benchmark points, the relic density, 
DM direct search, collider constraints and a strong FOPT 


4 345 can be very small due to the cancellation between three 
couplings A3, A4 and A5 while keeping A3 large enough to produce 
a strong FOPT. 


hz =1.68% at 240 GeV CEPC 
=0.0107 
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FIG. 1: The phase transition GW spectra h?Qaw in the IDM. 
The colored regions represent the expected sensitivities of 
GW interferometers U-DECIGO, DECIGO, BBO and eLISA, 
respectively. The black line depicts the GW spectra in the 
IDM for the set of benchmark points, which also represents 
the corresponding hZ cross section deviation at the 240 GeV 
CEPC and the corresponding DM coupling. 


can be satisfied simultaneously. Using the methods and 
formulae in the appendix, the phase transition GW signal 
is shown in Fig [i] which is just within the sensitivity of 
BBO and U-DICIGO. The colored regions represent the 
sensitivities of different GW experiments (DECIGO [53], 
eLISA [54], BBO, and U-DECIGO [B6]), and the black 
line corresponds to the GW signals, which also means 
the hZ cross section (et + e~ + h + Z) deviation from 
the SM in 240 GeV circular electron-positron collider 
(CEPC). At the 240 GeV CEPC with an integrated 
luminosity of 10 ab~', the precision of onz could be about 
0.4% [56]. And at the 240 GeV CEPC, the deviation of 
the hZ cross section 6,7 = ie at one-loop level [57] is 
about 1.68% [59], which is well within the sensitivity 
of CEPC. The international linear collider (ILC) can 
also test this model. The GW signal and the hZ cross 
section deviation at future lepton collider can make a 
double test on the DM of IDM as shown in Fig. 


III. GRAVITATIONAL WAVES IN MIXED 
DARK MATTERS 


We have discussed FOPT and GWs if there is only 
one single multiplet scalar dark matter in the dark sector. 
Due to the tight correlation between strong FOPT and 
the DM direct detection, only the IDM is viable for strong 
FOPT and detectable GWs. Based on the relic density 
requirement, the IDM has very limited viable parameter 
region: My ~ 55 ~ 75 GeV, and the blind spot region 
A345 ~ 0. In this section, we would like to extend the 
single DM multiplet models into the mixed scalar DM 
models. 

The mixed scalar DM scenario involves in several Z2- 


odd scalar multiplets in the dark sector, which could be 
mixed. The simplest models involve in two dark matter 
multiplets: the mixed singlet-doublet model (MSDM) 
and the mixed singlet-triplet model (MSTM) [89] [40]. 
Compared to the single scalar DM, the advantages of the 
mixed DM scenario are 


e It is easy to obtain a large blind spot region, at 
which the DM-Higgs coupling is zero.; 


e Ifthe DM is mixture of a singlet and multiplet (well- 
tempered), the relic density could be realized in a 
large DM mass range; 


e There are more degree of freedoms which contribute 
to the thermal barriers of the electroweak phase 
transition. 


From the above, we note that it is necessary to have a 
singlet component to reduce the large annihilation cross 
section during freeze-out. The higher multiplet n > 3 
is also tightly constrained by the plasma screening and 
Higgs diphoton constraints. 

We will focus on the strong FOPT and GWs in the blind 
spot region with broad DM mass range in the MSDM 
and MSTM. Here, we adapt the effective Lagrangian in 
Ref. [40]. Denote the SM singlet as S and the neutral 
component of the multiplet as H}. Due to the mixing 
between two fields, the mass eigenstates have 


x = coss —sindH?, (10) 
s = sins + cosOH?. (11) 

The effective Lagrangian for the mass eigenstates has 
Lot = Ghyyhxx + ansshss +++. (12) 


Define the dimensionless coupling Anyy = Gnyy/(2v) and 
Anss = 4nss/(2v). According to Ref. [40], the mixed DM 
scenario could have large viable parameter region after 
imposing the direct detection constraints and requirement 
on the DM relic abundance. Unlike the IDM, the DM 
mass could be arbitrary. The viable parameter region 
strongly depends on two parameters: the Anyy, and the 
mixing angle 0. The spin-independent direct detection 
rate is dominated by the Higgs exchange with the cross 
section: 


a2 m2, \? 
oe = FRB ( TN) (13) 
T Mymp 

Thus, the Aj. needs to be small to avoid tight direct 
detection constraints. In this study, we take the blind spot 
region: Anyy ~ 0. The dominant thermal annihilation 
channels could be very different and distinct depending 
on the DM mass region. To produce detectable GW 
signature, we would like to focus on the moderate DM 
mass region (mw < my < 300 GeV) 5. The dominant 
annihilation channels are 


5 Of course, the parameter space near the Higgs funnel region is 


e the s-channel annihilations:s yy > h > 
WW, ZZ, f f, in which the annihilation rate depends 
on the Any; 


e the contact four point annihilation: xx —> WW for 
real multiplet and yy —> WW, ZZ for doublet; 


e the ¢t-channel annihilation: yy > WW via exchang- 
ing of the charged scalars; 


Given Anyy ~ 0 in the blind spot region, we assume that 
the contact and t-channel annihilation processes are dom- 
inant °. In this case, the annihilation cross sections only 
depend on the gauge interactions, and the annihilation 
rate near the blind spot region is greatly simplified as 


> 3q* + 69? g? + g”? 


ae 
Oyy Vvu ~ 256mm? sin’ 0, (14) 
for the mixed singlet-doublet, and 
2 g“ 2 
Oxxswwv S 7s sin 0, (15) 
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for the mixed singlet-triplet model. To obtain the required 
relic abundance, the mixing angle sin 0 needs to satisfy 


Mx 


ing ~ — x _ l 1 

sin 0 540 CoV (doublet), (16) 
mMm 

inf ~ ——* — (triplet 17 

an 5000 Gay Se ey) 


which could be deduced from Ref. [37] [38]. Note that this 
simple relation is only valid near the blind spot region. 
In the interested mass region (mw < my < 300 GeV), 
the well-tempered DM is favored. Therefore, we could 
choose the benchmark points by acquiring Anyy is close 
to zero (blind spots) and sin @ is determined by the above 
Eq. (17). Using these two conditions, we could obtain 
correct relic abundance, while escaping the DM direct 
detection constraints and obtaining a strong FOPT with 
detectable GWs. 


A. The Mixed Singlet-Doublet Model 


The MSDM contains a singlet S with hypercharge Y = 
0, and a doublet with hypercharge Y = 1/2. The tree 
level potential in the MSDM is 


1 j 
v% = 5 Ms" + M2 H} Ho + 58S? ||? + A361 OH! Ap 


+ a|®t H|? + A (sim) +he]+A [SoH J he 


also viable. But we focus on a broader parameter region in this 
work. 

6 There is still small parameter region in which the s-channel and 
gauge interactions have interference effects. 


“esp © esm +2 (48)? + (22ta) 4 ny 


Here, we omit all interactions involving only S and Hə 
and the kinetic terms. The mixed mass matrices are 


T2 
2 Mg Av 
Neutral = ( 


Av M3 ) ? (18) 


where M2 = M2+ iv’ As and M2, = M3 + $V" Agas. The 
two eigenvalues of the mass-squared matrix are 


= = es = 2 
M + M2 \ (18 irg) +4v2A2| , (19) 


with the smaller eigenvalue corresponding to DM mass 
square 


7 k ` -42 
my = z | Ms + Mp — y (i12 - ig) + 4y?.A?| . 
(20) 
The mixing angle is 
2A 
tan 26 = ———__. (21) 
M2 — M? 
The coupling between the Higgs boson and DM is 
1 
anxx = 3Y (As + Asas) (22) 


2v A? + ty (m3 = MB) (As ami A345) 
F BEE 
y (13 - MB) +4024? 


and the coupling between the Higgs boson and the singlet 
scalar is 


1 
Ahaa = zr Os + 345) (23) 


2vA? + Lv (M3 — MB) (As — Asus) 
y(i - MB) + 42a? | 


In the DM favored parameter spaces (small Aryx, 
sind ~ m,/540 GeV and mw < my < 300 GeV), a 
strong FOPT can be induced due to the fact that the 
strong FOPT can occur in the IDM and the mixing with 
an inert singlet can enhance the FOPT. The mixing term 
A [S@Di + h.c.] does not significantly change the phase 
transition property except for the change of the eigenval- 
ues of the masses. Taking the same approximations as 
in the IDM and writing in the forms of Eq. (5), we have 
ye 4 
E c. The strong FOPT can be easily realized if A3 
and Aass are order one which are allowed by the above 
DM constraints. 

We take one set of benchmark points A3 = 3.006, Ay = 
—1.5, As = —1.5, Ag = 4.006, A = 91 GeV and Mp = 
Ms = 117.5 GeV. Then, we have mẹ = 100.6 GeV, 


+ 
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FIG. 2: The phase transition GW spectra h?Qeaw in the 
MSDM. The colored regions represent the expected sensitivi- 
ties of GW interferometers U-DECIGO, DECIGO, BBO and 
eLISA, respectively. The black line depicts the GW spectra 
in the MSDM and the corresponding DM coupling for the set 
of benchmark points. 


Aryx = 9.0009 and Anss = 2.005. The corresponding GW 
spectra are shown in Fig [2] This set of benchmark points 
can evade the constraints from DM direct experiments, 
give the correct relic density and produce detectable GW 
signals by eLISA, BBO, DECIGO and U-DECIGO. We 
estimate the deviation of the hZ cross section 6;,z at one- 
loop level is about 2.78% [58] [59] at the 240 GeV CEPC, 
which is well within the sensitivity of CEPC and ILC. As 
for the precise prediction on the collider signals at the 
CEPC and ILC in this MSDM model, we leave it in our 
future study. 


B. The Mixed Singlet-Triplet Model 


Similarly, the MSTM contains a real singlet and a real 
triplet with Y = 0. The relevant potential in the MSTM 
is 

1 

V = 5 Mss + M3Tr(H3) + ky! ®Tr( HZ) 
n “japs + ESD’ H,®. (24) 

The coupling between the Higgs boson and DM is 

1 
anx = 3” (k+ Ky) (25) 
tv3? + gu Ce Wig) (K — Ky) 
7 zayi 
y (123 - M) + pote? 


with the DM mass square 


? 


a 
2 aZ 
My = 5 


A , (26) 
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FIG. 3: The phase transition GW spectra h?Qew in the 
MSTM. The colored regions represent the expected sensitivi- 
ties of GW interferometers U-DECIGO, DECIGO, BBO and 
eLISA, respectively. The black line depicts the GW spectra 
and the corresponding DM coupling in the MSTM for the set 
of benchmark points. 


where M2 = M2 + iv? and MZ = M2 + *2v?. The 


coupling between the Higgs boson and the singlet is 


1 
Ghss = 5¥ («+ Ky) (27) 


pus? + sv (2 Wiz) (K — Ky) 
-mre 


Taking the same approximations as in the above discus- 
sions, we have eusrm = ESM + Ord + (Anss)? + 
2 (ky /2)?/. In the MSTM, we take the set of bench- 
mark points « = 0.01, ks = 3.0, € = 0.31337, My = 50 
GeV, Ms = 119.93 GeV. Then, we have m, = 120 GeV, 
Anxx = 9.001 and Anss = 1.504. The GW spectrum is 
shown in Fig[3] which is within the sensitivity of the 
BBO, eLISA, DECIGO and U-DECIGO. We also leave 
the study of collider signals in the MSTM in our future 
study. 


IV. CONCLUSION 


We have investigated the GW signatures associated 
with the possible strong FOPT in the early universe in 
a generic class of inert scalar multiplet DM models. We 
have found that there is usually a strong tension between 
the strength of the phase transition and the DM direct 
detection, except the blind spot region with narrow DM 
mass range in the IDM. In the IDM, we have shown that 
although the direct detection experiments could not probe 
the blind spot region, at which the Higgs-DM coupling is 
zero, the GW detection and lepton colliders could help 


us further explore the parameter space which cannot be 
addressed by the DM searches. 

To enlarge the blind spot region with broader DM 
mass range, we considered the mixed scalar DM scenario. 
We focused on the minimal singlet-doublet dark matter 
and minimal singlet-triplet DM models, and showed that 
the blind spot region in these models could avoid all 
the relevant constraints, have required relic density, and 
produce strong FOPT in the early universe. The GW 
signatures are explored in this scenario, and we found 
that the blind spot region could be further explored via 
the GW signatures and collider signatures in future. This 
provides us a new way to detect the DM in future. And our 
strategy could also be used to more general DM scenarios 
than the scalar multiplet DM models. We expect the 
future GW and collider experiments, such as eLISA and 
TianQin [60], CEPC and ILC, could probe these DM 
scenarios beyond the reaches of the DM direct detection 
experiments. 
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Appendix 


In this appendix, the details on the calculation of phase 
transition GW spectrum are shown. In general, there ex- 
ist three mechanisms to produced phase transition GWs: 
the first one is the well-known bubble collisions [5], the 
second one is the turbulence in the fluid, where a certain 
fraction of the bubble walls energy is converted into tur- 
bulence [7] [8], and the last one is the new mechanism of 
sound waves [6]. 

When a strong FOPT occurs, via thermally fluctuating 
or quantum tunneling the potential barrier, bubbles are 
nucleated with the nucleation rate per unit volume [ = 
To(T)e~ 8") and To(T) « T* [61]. The Euclidean action 
Se(T) ~ S3(T)/T [62] [63], and then T = Pye 98/7 


where 
S(T) = f ae [509 + Vale, T) . (28) 


Thus, the bubble nucleation rate can be obtained [64] by 
calculating the profile of the scalar field through solving 
the following bounce equation: 


PË 2d® IVg(®,T) 


T — 2 
dr? rdr O® 9, (a 


with the boundary conditions 4? (r = 0) = 0 and ®(r = 
oo) = 0. The phase transition terminates if it satisfies 


S3 (T) /T; = 41n(T, /100GeV) + 137. (30) 


The phase transition GW spectrum depends on four pa- 
rameters. The first parameter is a = es. where T, is 
determined by Eq. (80). The false vacuum energy (latent 
heat) density e(T,) = [TZE — væn (T)]|r=r, , and the 
plasma thermal energy density praa(Tx) = Tg, Ga igs 
The parameter a represents the strength the FOPT, 
namely, a larger value of œ produces stronger GW sig- 
nature. The second parameter is # = TLS) ; 
: T=T, 
where one has 6 = -Be ji, ~ idli, namely, 67! 
corresponds to the typical time scale of the phase tran- 
sition. The third parameter is the efficiency factor A; 
(i=co,tu,sw), and the last parameter is the bubble wall 
velocity up. 

Once the four parameters are known, the corresponding 
phase transition GWs in the three mechanisms can be 


directly obtained after considering the red-shift effects 


t 


1/6 
ax = 1.65 x 107°Hz x T ( T, \(5) . The current 


ao ToOGev 

peak frequency at each mechanism is fi = f/a./ao with 
i=co,tu,sw, respectively. In the bubble collision mecha- 
nism, the GW spectrum is expressed as [65] 


Qeo( fh? ~1.67 x 10-9(42) "(Aaa)" (2)! 


Iiv : (f/ feo)?’ 
ü (ap sas) f J] 


with the peak frequency f*, = 0.628/(1.8—0.1v,+v2) [65] 
at Ta. In the turbulence mechanism, the GW signals have 
the peak frequency at about ff, = 1.758/w at T, [54], 
and the phase transition GW intensity is formulated by 


[8] [66] 
Wu (f)h? =3.35 x 10-4(=) (wey) a 


: (Ffr)? 
OF Ffa) + 8rfao/0 H) 


In the sound wave mechanism, the GW spectrum can be 
written as [6] 


Oy (A)? 2.65 x 10-8 (He) (S2)? (100) 


B/\1+a gi 
TF /Fow)9!7 77/2 
i EEA. 


with the peak frequency fă, = 28/(V3u,) at T. [6l BA, 
where A,,, ~ a (0.73 + 0.083,/a + a)’ for relativistic 
bubbles [67]. 
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